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Abstract. Pure N-body high-resolution re-simulations of 13 massive dark mat- 
ter halos in a ACDM cosmology are presented. The resulting sample is used to 
investigate the physical origin of the Fundamental Plane, "Faber- Jackson" , and 
"Kormendy" relations observed for nearby galaxy clusters. In particular, we focus 
on the role of dissipationless hierarchical merging in establishing or modifying 
these relations. Contrarily to what found in the case of the Faber-Jackson and 
Kormendy relations for galaxies (see Londrillo, Nipoti & Ciotti, this conference), 
dissipationless merging on cluster scales produces scaling relations remarkably 
similar to the observed ones. This suggests that gas dissipation plays a minor role 
in the formation of galaxy clusters, and the hierarchical merger scenario is not at 
odd with the observed regularity of these systems. 
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1. Introduction 

Like early- type galaxies, also galaxy clus- 
ters follow well defined scaling relations 
among their main observables: in partic- 
ular, a luminosity-velocity dispersion re- 
lation (i.e., a "Faber- Jackson-like" rela- 
tion), a luminosity-radius relation (i.e., 
a "Kormendy-like" relation), and the 
Fundamental Plane (hereafter FP; e.g., 
Schaeffer et al. 1993, hereafter S93; Adami 
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et al. 1998; and for ellipticals: Djorgovski 
et al. 1987; Dressier et al. 1987; Faber 
& Jackson 1976; Kormendy 1977). Besides 
their potential importance as distance indi- 
cators, these relations contain information 
on the cluster formation and evolution pro- 
cesses. Within the commonly accepted cos- 
mological scenario, where cold dark matter 
(CDM) is the dominant component of the 
Universe, structure formation is driven by 
hierarchical dissipationless merging: small 
DM halos form first, then they merge to- 
gether to form larger systems, and the evo- 
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lution of the baryonic matter follows that 
of the hosting DM halos. 

The importance of investigating the 
physical origin of the observed scaling re- 
lations, and the role of dissipationless hi- 
erarchical merging in establishing or mod- 
ifying them, is therefore apparent, and ex- 
plains the large number of theoretical works 
devoted to this problem in the case of el- 
liptical galaxies (e.g. Capelato et al. 1995; 
Nipoti, Londrillo, Ciotti 2002a, b; Gonzales 
& van Albada 2002; Evstigneeva et al. 
2002; Dantas et al. 2002; Londrillo et al., 
this conference). In particular, it has been 
shown that, while the FP is reproduced 
by the end-products of equal mass merg- 
ers, the Faber-Jackson and the Kormendy 
relations are not (Nipoti et al. 2002ab; 
Londrillo et al., this conference). In the 
case of galaxy clusters however, almost 
no theoretical studies are devoted to their 
FP (see, e.g., Pentericci, Ciotti & Renzini 
1995; Fujita & Takahara 1999; Beisbart, 
Valdarnini & Buchcrt 2001). 

We address this problem in the present 
work, by means of high resolution re- 
simulations of very massive DM halos (the 
hosts of galaxy clusters) in a ACDM cos- 
mology. Following their hierarchical dis- 
sipationless evolution in detail, we verify 
whether they reproduce the observed scal- 
ing relations, both at the present time, and 
at higher redshift. A full analysis of this 
work will be given elsewhere (Lanzoni et 
al., in preparation), while we present here a 
selection of the main astrophysical results, 
focusing on the numerical issues. 

2. High-resolution re-simulations 

To investigate whether the dark mat- 
ter counterpart of galaxy clusters, as ob- 
tained by numerical simulations, do de- 
fine the observed scaling relations, a large 
enough sample of very massive DM ha- 
los was needed. For this purpose, we have 
employed the dissipationless "Very Large 
Simulations" (VLS; see Yoshida, Sheth, & 
Diafcrio 2001), where the simulated comov- 
ing volume of the Universe is sufficiently 



large: the box side is of 479 h 1 Mpc, with 
H = 100 h,- 1 km s" 1 Mpc" 1 , and h = 
0.7. The adopted cosmological model is a 
ACDM Universe with fl m = 0.3, fl A = 0.7, 
spectral shape T = 0.21, and normalization 
to the cluster local abundance, er 8 = 0.9. 
The total number of particles is 512 3 , of 
6.86 x \{) w M Q /h mass each. 

From these simulations, we have se- 
lected a sample of 13 halos with masses 
between 10 14 M Q //i and 2.3 x 10 15 M Q /h. 
They span a variety of shapes, from nearly 
round to more elongated. The richness of 
their environment also changes from case 
to case, with the less isolated halos usu- 
ally surrounded by pronounced filamen- 
tary structures, containing massive neigh- 
bors (up to 20% of their mass). All the se- 
lected halos are required to have a mas- 
sive progenitor at redshift z = 0.5 (their 
masses range between 6 x 10 13 and 1.5 x 
10 15 M Q //i), and also at z = 1 (with masses 
between 4 x 10 13 and 7 x 10 14 M Q //i). Such 
a choice allows us to study galaxy clusters 
not only locally, but also at higher redshift. 

Given the mass resolution of the sim- 
ulations, less than 1500 particles compose 
a halo of 10 14 M Q //i, and its properties en- 
tering the FP relation cannot be accurately 
determined. We have therefore resimulated 
at higher resolution the halos in our sam- 
ple, by means of the technique introduced 
by Tormcn, Bouchct & White (1997), that 
we briefly summarize in the following. 

The first step is to select, in a given cos- 
mological simulation, the halo one wants 
to reproduce at higher resolution. Then, 
all the particles composing the selected 
halo and its immediate surroundings are 
detected in the initial conditions of the 
parent simulation: the number of particles 
within the region thus defined is increased, 
and the region is therefore called "the 
high resolution (HR) region". As a con- 
sequence, the mean inter-particle separa- 
tion decreases, and the corresponding high- 
frequency modes of the primordial fluctua- 
tion spectrum are added to those on larger 
scales originally used in the parent sim- 
ulation. The overall displacement field is 
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also modified consequently. At the same 
time, the distribution of surrounding par- 
ticles is smoothed by means of a spherical 
grid, whose spacing increases with the dis- 
tance from the center: in such a way, the 
original surrounding particles are replaced 
with a smaller number of macroparticles, 
whose mass grows with the distance from 
the HR region. Thank to this method, even 
if the number of particles in the HR re- 
gion is increased, the total number of par- 
ticles to be evolved in the simulation re- 
mains small enough to require reasonable 
computational costs, while the tidal field 
that the overall particle distribution exerts 
on the HR region remains very close to the 
original one. For the new initial configura- 
tion thus produced, vacuum boundary con- 
ditions are adopted, i.e., we assume a van- 
ishing density fluctuation field outside the 
spherical distribution of particles with di- 
ameter equal to the original box size L. A 
new N-body simulation is then run start- 
ing from these new initial conditions, and 
allows to re-obtain the selected halo at the 
suited resolution. 

We have applied this technique to the 
13 massive DM halos selected in the VLS. 
For 8 out of the 13 selected halos, the 
resolution has been increased by a fac- 
tor ~ 33, using high-resolution particles 
of - 2.07 x l0 9 M Q /h each, while a fac- 
tor of 2 improvement has been adopted for 
the 5 intermediate mass halos (the parti- 
cle mass is 10 9 Af Q //i in this case). The 
gravitational softening used for the high- 
resolution region is e = 5kpc//i, corre- 
sponding to about 0.2% and 0.5% of the 
virial radius of most and less massive ha- 
los, respectively. This scale length repre- 
sents the spatial resolution of the resimula- 
tions, to be compared with that of 30 kpc/h 
of the original VLS. Note that to prevent 
low-resolution macroparticles to "contam- 
inate" the final resimulated halo (i.e., to 
end within its virial radius), we need to de- 
fine the HR region not only through the 
particles of the halo itself in the parent 
simulation, but considering also a bound- 
ary region in its immediate surroundings. 



This is particularly important in the case 
of less massive halos, since the small num- 
ber of their composing particles does not 
allow to define their encompassing region 
in the initial conditions with enough pre- 
cision. To avoid such a "contamination" 
problem in our re-simulations, a radius of 
about 3 r v i r and 5 r v j r for the boundary re- 
gion was necessary for the most and the less 
massive halos, respectively. The resulting 
number of high and low resolution particles 
are listed in Table 0, together with their 
sumQ. To run the resimulations, the par- 
allel dissipationless tree-code "GADGET" 
(Springel, Yoshida & White 2001) has been 
used on the IBM SP2 and SP3 of the Centre 
Informatique National de PEnseignement 
Superieur (CINES, Montpellier, France), 
and the CRAY T3E of the RZG Computing 
Center (Munich, Germany) , that have com- 
parable processor speed. The number of 
processors and the CPU time required for 
the corresponding runs are also listed in 
Table [J 

Results have been dumped in output 
files for 100 time steps, equally spaced in 
the logarithm of the expansion factor of 
the Universe between z = 20 and 2 = 0. 
At each time step, a Spherical Overdensity 
(Lacey & Cole 1994) halo finder has been 
used to detect the halos, and to estimate 
their virial mass and radius (see Table |l|). 
With respect to the original virial masses 
and radii, those of the resimulated halos 
show non-systematic and small differences 
(of the order of few percent), particularly 
for the most massive objects. Moreover, 
also the halo formation history (the way 
how mass assembly proceeds with time) is 
very similar to that in the parent simula- 
tion, thus assuring the reliability and pre- 
cision of the high-resolution re-simulation 
technique we use. Visually, the noticeable 
improvement provided by such a technique 
is apparent from Figs. |l| and ^, where the 



1 Note that to get the same mass resolution 
(2.07 x W 9 M Q /h or 10 9 M Q //i) in the entire 
volume of the parent simulation, about 1650 3 
of 2100 3 particles of would have been required. 
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images of the 3 most massive halos, before 
and after the resimulations, are shown at 
the scales of ~ 30Mpc/ft. and ~ 2r V j r , re- 
spectively. 

3. Observed scaling relations at z = 

Analyzing a sample of 16 nearby galaxy 
clusters, S93 found the following relation 
between the observed luminosity and ve- 
locity dispersion: 

L oc a 1 ' 8 ™ 4 , (1) 

where L is in units of L* = 1.325 x 
10 10 /i _2 L© (h — 1), and a in km/s, and 
between luminosity and radius: 

L oc i?i-34±o.u (2) 

where R is the effective radius measured 
in Mpc. In both cases, the observed scat- 
ter is very large, while a considerable im- 
provement is found when linking the three 
observables together in a FP-like relation. 
While a bi-paramctric fit was performed by 
S93, here we derived the FP by means of 
the Principal Component Analysis (here- 
after PC A; see e.g. Murtagh & Heck 1987), 
that searches for the most suitable combi- 
nations of the three observables able to de- 
scribe a thin plane. By performing a PCA 
on the data sample of S93, we obtain the 
new orthogonal variables 

Pi = a t log R + [3, log L + 7 l log a, (3) 

with R in Mpc/h, L in 100L*, a in 1000 
km/s. In Fig. |^ we show the resulting dis- 
tribution of the observed clusters (blue tri- 
angles) in the {pi,P3) and (pi,P2) spaces, 
the former providing an exact edge-on view 
of the FP and making apparent its small 
thickness. 

4. Scaling relations at z = for 
simulated clusters 

To derive the analogous quantities for the 
simulated clusters, we have constructed 
their projected radial profiles by counting 
the DM particles within concentric shells 



around the center of mass, for three orthog- 
onal directions. Thus, the projected half- 
mass radius Rh has been derived as the ra- 
dius of the shell containing half the total 
number of particles, and the velocity dis- 
persion a has been computed by averaging 
the squared line of sight component of the 
(barycentric) velocity over all the particles 
within Rh. Since the simulated clusters, as 
well as the real ones, are not spherical, such 
a procedure gives different values of i?h and 
a for the three considered line of sights, 
the maximum difference never exceeding 
33% and 21% in the two cases, respectively. 
Finally, for converting the DM mass into 
light, we have adopted a mass-to-light ra- 
tio with a small dependence on the lumi- 
nosity, as suggested by observations (S93; 
Girardi et al. 2002): M/L oc L a , with 
a = 0.3. In particular, we have assumed: 
M/L = 450(M/Af ) Q /( Q+1 )M Q /L Q , where 
M = 7.5 x 10 14 Af Q . 

By performing a PCA on the 13 DM ha- 
los for the three considered lines of sights, 
using i?h, c, and L computed as just de- 
scribed, we obtain a very well defined and 
thin FP. To check whether this FP is sim- 
ilar to the observed one, we have also con- 
structed the PCA variables pi by combining 
i?h, c, and L derived for the DM halos, with 
the coefficients ai, Pi, ji obtained for the 
observed clusters. The resulting FP is prac- 
tically indistinguishable from the observed 
one, as apparent from Fig.|^. Moreover, also 
the L-a and L-R relations of simulated and 
real clusters are in remarkable agreement 
(Fig. |). 

Note that assuming a constant M/L ra- 
tio, the FP of simulated clusters is still well 
defined, and only shows a mild tilt with re- 
spect to the observed one (see solid line in 
Fig. |). 

5. Scaling relations at higher redshift 
for simulated clusters 

While no observational data are yet avail- 
able for galaxy clusters at higher redshift, 
our high-resolution re-simulations allow to 
investigate the scaling relations of their 
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Table 1. Characteristics of the high-resolution resimulations of the 13 massive halos: 
Nhr = number of high-resolution particles, Nlr = number of low-resolution macropar- 
ticles, Nt = total number of particles, N proc = number of processors, icpu = hours per 
processor required for the resimulations, M v i r = virial mass in 10 14 M Q //i, r v i r = virial 
radius in Mpc/h. 



Name 


N HR 


N LR 


N T 


Nproc 


icpu 


M vir 


/"vir 


g8 


3700120 


191733 


3891853 


32 


80 


23.42 


2.75 


gl 


2574717 


202301 


2777018 


32 


70 


13.99 


2.31 


g72 


3299865 


194277 


3494142 


32 


57 


11.77 


2.18 


g696 


4870197 


184314 


5054511 


64 


60 


11.37 


2.16 


g51 


1677364 


213477 


1890841 


32 


44 


10.78 


2.12 


g245 


3437317 


215269 


3652586 


16 


18.4 


6.50 


1.79 


g689 


3252085 


215809 


3467894 


16 


18.3 


6.08 


1.75 


g564 


2068981 


227187 


2296168 


16 


12.8 


4.91 


1.63 


gl777 


3094123 


219047 


3313170 


16 


18.0 


3.83 


1.50 


g4478 


2293433 


225706 


2519139 


16 


12.8 


2.92 


1.37 


g914 


250605 


247091 


497696 


16 


7.3 


1.45 


1.09 


g3344 


206140 


248756 


454896 


16 


6.3 


1.09 


0.99 


g!542 


207202 


248948 


456150 


16 


6.3 


0.82 


0.99 



DM counterparts at any epoch. For that 
purpose, the most massive progenitors at 
z = 0.5 and at z = 1 of our 13 DM ha- 
los have been considered, and their pro- 
jected half-mass radius, velocity dispersion 
and luminosity have been derived with the 
same procedure used for the z = objects. 
Also in this case, when combining i?h, cr, 
and L with the coefficients on, Pi, 7^ de- 
rived for the real clusters at z = 0, the re- 
sulting FP is very similar to the edge-on FP 
observed locally, while when seen face-on, 
the region populated by the high-z DM ha- 
los is systematically shifted towards larger 
(smaller) values of p\ (j>2) with respect to 
the region occupied by the observed nearby 
clusters. In addition, very well defined L- 
a and L-R relations are already in place 
at these redshifts. With respect to the ob- 
served L-a relation at z — 0, a flatter slope 
is found for increasing rcdshift for the sim- 
ulated clusters, while the opposite trend is 
found in the case of the L-R relation, that 
becomes steeper for increasing z. 



6. Discussion and conclusions 

We have described a technique for increas- 
ing the (mass and spatial) resolution of ob- 
jects selected in a given cosmological sim- 
ulation. Applying it to sample of 13 mas- 
sive halos at z = in a ACDM cosmology, 
we have investigated whether the scaling 
relations observed for nearby galaxy clus- 
ters are also followed by their DM hosts, in 
a dissipationless hierarchical merging sce- 
nario. By considering the most massive pro- 
genitors at z = 0.5 and z = 1 of the 13 
selected halos, we have also investigated 
whether the same scaling relations were al- 
ready in place at those earlier epochs. 

The main conclusions can be summa- 
rized as follows: 

— the high-resolution re-simulation tech- 
nique is very powerful and reliable. The 
precision of the method depends on the 
number of particles composing the orig- 
inal object, and on the dimension of the 
boundary region considered to avoid the 
"contamination" from macroparticles. 
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— The DM hosts of galaxy clusters do 
define a FP that is practically indis- 
tinguishable from the one observed for 
nearby galaxy clusters. 

— A very good agreement between the 
scaling relations of simulated and ob- 
served clusters is also found in terms of 
the L-a and the L-R relations. 

— The FP of simulated clusters is already 
defined at z — 0.5 and z = 1, and it is 
very similar and as thin as the one ob- 
served locally; the only differences are 
due to an obvious shift towards smaller 
radii and masses for increasing z. Also 
the L-a and L-R relations are already 
in place at high redshifts, with a differ- 
ence in the slope (with respect to the re- 
lations observed locally) due to the fact 
that at recent epochs, more massive ob- 
jects evolve more rapidly than smaller 
systems. 

All the results about the scaling re- 
lations of simulated and observed clus- 
ters are obtained by assuming a mass-to- 
light ratio with mean value of the order of 
450hM Q I Lq (well in agreement with what 
estimated from observations of galaxy clus- 
ters), and with a dependence on the lumi- 
nosity also suggested by the observations. 

Given the very different nature of the 
two components (baryonic and non bary- 
onic) and of the physical processes acting 
on them, the fact that the DM halos define 
the same relations shown by real clusters is 
not obvious at all. Moreover, dissipation- 
less hierarchical merging at galactic scales 
is unable to reproduce the observed scal- 
ing relations of elliptical galaxies: the end- 
products of major mergers do lie on the ob- 
served FP, but substantial deviations from 
the Faber-Jackson and the Kormendy rela- 
tions are found in this case (see Londrillo et 
al., this conference). Such a result suggests 
that, while dissipation has a major role in 
the formation and evolution of ellipticals, it 
is negligible with respect to the pure grav- 
ity in settling the properties of galaxy clus- 
ters. 
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Fig. 1. Images of the 3 most massive DM halos in the parent simulation (left panels) 
and after the high-resolution re-simulation (right panels). All panels show the projection 
along an arbitrary line of sight, of a region of ~ 30Mpc//i around the halo center of 
mass. 




Fig. 2. Images of the 3 most massive DM halos in the parent simulation (left panels) 
and after the high-resolution re-simulation (right panels). All panels show the projection 
along an arbitrary line of sight, of a region of rN -' Z Tvir 

around the halo center of mass. 
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Fig. 3. Observed (blue triangles) and simulated (red hexagons) clusters in the edge-on 
and face-on view of the FP (upper and lower panels, respectively), at redshift z = 0. The 
linear best fit to the simulated clusters in case of a constant M/L is shown by the pink 
solid line. 
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Fig. 4. L-a (left panels) and L-R (right panels) relations for the observed (blue triangles 
and dashed line) and simulated (red hexagons and solid line) clusters, at redshift z = 0. 
Luminosities are in units of 100 L*, velocity dispersions in km/s, and radii in Mpc. 



